Structure of and E2 transition in 16 C in a 14 C+n+n model 
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Abstract 

A three-body model of 14 C+n+n is applied to study the energy spectrum and the hindered E2 
transition in 16 C. A realistic two-nucleon potential is used for the valence neutrons. Both spin 
singlet and triplet components for the neutrons are taken into account. The three-body problem 
with a Pauli constraint is solved in a stochastic variational method. For the n- 14 C potential chosen 
to reproduce the properties of 15 C, the low-lying energy spectrum agrees reasonably well with 
experiment, but the ground state is predicted to be about 1 MeV high. The calculated B{E2; 
2^ ^O] 1 ") value is about twice the measured value if the polarization charge of the valence neutrons 
is taken to be the same as that required to fit the 15 C data. The correlated motion of the valence 
neutrons is displayed through the two-neutron density distribution. 
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Recently the E2 transition from the first 2 + state to the ground + state in 16 C has been 
studied through a lifetime measurement using a recoil shadow method 1] and 16 C+ 208 Pb 
inelastic scattering Q. The B(E2) value is found to be 0.63 ±0.12 e 2 fm 4 , which corresponds 
to anomalously small strength of about 0.26 W.u. The anomaly is apparent by comparison 
with other C isotopes. The B(E2; 2± — > 0^) values of 10 C, 12 C and 14 C decrease gradually as 
expected from the increasing excitation energies. However, the case of 16 C strongly deviates 
from this expectation: The excitation energy of the 2f state of 16 C is only 1.77 MeV, much 
smaller than the case of 14 C (7.01 MeV), but the B(E2) value of 16 C is nevertheless much 
smaller than that of 14 C, which is 3.74 e 2 fm 4 (1.87 W.u.). The B(E2) value for the -> 
transition in 15 C is also small, 0.97 ± 0.02 e 2 fm 4 (0.44 W.u.) 0. 

In the previous paper jij], a three-body model of 14 C+n+n was applied to study the 
anomalous E2 transition in 16 C under the assumption that the 14 C core has a neutron shell 
closure of the p shell and the two valence neutrons are in a spin-singlet state. The 14 C+(sd) 2 
model for the low-lying levels of 16 C is supported by experiment as well jHEf. The calculated 
B(E2) value was 1.38 e 2 fm 4 , 2.2 times larger than experiment if one uses the same neutron 
polarization charge as that required for 15 C. The deformations of both proton and neutron 
densities of the C isotopes are studied in antisymmetrized molecular dynamics calculation [?J , 
where the neutron shape is found to change significantly with the neutron number and the 
B(E2) value of 16 C in a variation after projection calculation is 4-6 times too large compared 
to experiment. The B(E2) value is overestimated also in the deformed Hartree-Fock model 
and the shell model calculations (If. The purpose of the present investigation is to study 
the structure of 16 C thoroughly and to reexamine a mechanism which leads to the hindered 
transition. Our basic assumption is the same as before, so the relevant levels of 16 C are 
assumed to be generated from the 14 C+n+n model. There are, however, some noticeable 
differences. We here use more general basis functions allowing for both spin-singlet and 
spin-triplet neutrons. This makes it possible to obtain unnatural parity states such as 3 + . 
In addition, we use a realistic interaction for the two neutrons in order to avoid uncertainties 
of the model, while in the previous paper an effective interaction was used and its strength 
was increased so as to fit the ground-state energy. 

The wave function for 16 C is determined from the following Hamiltonian 

H = T R + T r + U 1 + U 2 + v 12 , (1) 

where the subscripts of the kinetic energies stand for the relative distance vector, R, from 
the center-of-mass of C to that of the two neutrons and the relative distance vector of the 
neutrons, r. The two-neutron potential vrz is taken from G3RS (case 1) potential which 
contains central, tensor and spin-orbit forces and reproduces the nucleon-nucleon scattering 
data as well as the deuteron properties. The n- 14 C potential U takes the form 

U = -V f(r)+V 1 £-s~f(r), (2) 

r ar 

i 

where /(r) = [l + exp( Iz ^)]~ 1 with R c =r A(? (A c = 14). The parameters of U are determined 

to reproduce the energies of the L and | x states of 15 C. In order to make the lsi/2 state 
lower than the 0d 5 /2 state, the £■ s strength V\ was chosen in the previous study to be about 
half of the standard value jicj ]. This potential is called set A hereafter. We also test other 
possibilities to fit these single-particle (s.p.) energies. In set B we use the standard value of 
V\ but weaken Vq for all the partial waves other than s-wave, while in set C a potential with 
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TABLE I: Parameters of n- C potential. 





Fo(MeV) 


Fi(MeVfm 2 ) 


r (fm) 


a(fm) 


set A 


-50.31 


16.64 


1.25 


0.65 


set B 


-50.31 (J = 0), -47.18 (Z^O) 


31.25 


1.25 


0.65 


set C 


-51.71 


26.24 


1.20 


0.73 



a larger diffuseness parameter is chosen. The parameters of each set are listed in Table UJ 
These potentials generate almost the same s.p. wave function for both lsi/ 2 and 0d 5 / 2 . 
Trial wave functions are assumed as a combination of correlated Gaussian bases 



K 



*jm(1,2) = y ZC i V JM (\i,A i ), (3) 

i=i 

* JM (\,A) = (l-P 12 ){e-^ a! [[%(a ;i )^ 2 (a !2 )] L X5(l,2)]jM}, (4) 

where Pu is a permutation of the two neutrons. The basis function is specified by a set of 
angular momenta A = (£%, £ 2 , L, S) and a 2x2 positive-definite, symmetric matrix A. Here 
xAx stands for Auxl+2Ai2XyX2+ A22X2 , and X\=R+\r and X2 =R—\r are the distance 
vectors of the neutrons from the center-of-mass of the 14 C core. The two neutrons are 
explicitly correlated due to the presence of the cross term A^xi-x-i, the inclusion of which 



is very important to obtain a precise solution |ll|, 12]. The angular parts of the two-neutron 



orbital motion are described using D^ m (r) =r e Y (m (r) and they are coupled with the spin 
part xs to the total angular momentum J. In the previous study, the angular motion was 
described with much simpler functions specified by a global vector but here a general 
form is used to take into account important L, S contents of the wave functions. 

The energy and the corresponding wave function are determined from a solution of the 
generalized eigenvalue problem 

YXHy - EB i3 \C 3 =0 (i = 1, 2, . . . , K), (5) 
(VM^Ai) |( f )\^ JM (X j ,A j )). (6) 



Hij 



It is vital to take into account the Pauli principle for the motion of the valence neutrons. 
Under the assumption that 14 C is neutron-closed, the Pauli constraint is fulfilled by imposing 
that the trial wave function has no overlap with any orbit u n ij m occupied in the 14 C core 

( Untjm (i)\^j M (l,2)) = (^ = 1,2), (7) 

where the s.p. orbit u n tj m is generated from U and n£j runs over Osi/2, 0^3/2, and Opi/2. 
The condition (J7|) is practically achieved by the orthogonal projection method The 
accuracy of the present calculation is such that the probability of mixing-in of the occupied 
orbits is of the order of 10~ 4 . 

Before superposing basis functions with different A, we did a pilot calculation in a single 
A channel. We use the algorithm called the stochastic variational method (SVM) 3 
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to optimize the parameter matrices A. The SVM increases the basis dimension one by one 
by testing a number of candidates which are chosen randomly, and in addition fine-tunes 
the already chosen parameters by a refinement process. The basis selection with the SVM 
is effectively performed not only to take care of the short-range repulsion of v 12 but also 
to satisfy the condition (J7J). For J 7r =0 + , those channels which give an energy lower than 
the 14 C+n+n threshold are only A=(0, 0, 0, 0) and (2, 2, 0, 0), while for J 7T =2 + such a bound 
solution is obtained only for A=(0, 2, 2, 0). Channels with > 3 play a minor role. We also 
studied levels with J 7r =3 + and 4 + . The most important channel is A=(0,2,2, 1) for 3 + and 
A=(2, 2, 4, 0) and (2, 2, 3, 1) for 4 + , respectively. Thus the 5=1 channel is important for 4 + . 

The basis sets chosen in the single-channel calculation serve as the basis functions for 
a coupled-channels calculation. Channels are truncated to those with £\ + £2 <4. A basis 
dimension K used in the calculation is about 450 for + and 500-650 for the other cases. 
Figure n compares the calculated energy levels of 16 C with experiment. The [/-dependence 
of the spectrum is moderate. The ground-state energy is about 1 MeV too high compared 
to experiment, but the energies of the other J n states are in fair agreement with experiment. 
No state with J 7T =1 + was obtained below the 15 C+n threshold. 

Table |H] summarizes some properties of the states calculated with set B potential. The 
other potentials give similar results. The root mean square radius of 16 C, r ( 16 C) , for the point 
nucleon distribution is calculated from the equation, r 2 ( 16 C) = ^r 2 ( 14 C) + ^(i? 2 ) + 2§g(r 2 ), 
where we substitute the point proton radius, 2.35 fm [3J| for r( 14 C). The r( 16 C) value for the 
ground state agrees well with that extracted from the reaction cross section analysis [15]. 
The probability of spin-singlet components, Ps=o, shows that our previous calculation |4( 
with no S = 1 components is not very good for the 2f state but quite acceptable for the 
ground state. P ss , P s d and Pdd stand for the probabilities for the two neutrons to occupy the 
(lsi/2) 2 , (isi/ 2 0d 5 /2) and (Oc^) 2 configurations, respectively. The sum of the probabilities, 
P ss +P s d+Pdd, is less than unity, which signals the importance of unbound s.p. orbits or n+ 14 C 
continuum states. Its effect appears more important in the 0^ and 2f states. The previous 
calculation gave P ss = 0.49 and Pdd = 0.39 for the ground state, which is consistent with the 
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FIG. 1: Low-lying energy levels of C for different sets of the n- C potential. Energy is given 
from the 14 C+n+n threshold. 
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TABLE II: Properties of the states in C. Set B parameters are used. Energy and length are in 
MeV and fm, respectively. 



r 


Ecal 


E 


<*?> 


(R 2 ) <r 2 ) 


(x r x 2 ) 


r( 16 C) P s= o 


p 

± ss 


Psd Pdd 


of 


-4.40 


-5.469 


17.8 


10.4 29.5 


3.04 


2.63 


0.92 


0.52 


0.37 


0+ 


-2.52 


-2.466 


19.7 


10.4 37.3 


1.10 


2.67 


0.86 


0.43 


0.53 




-3.18 


-3.699 


16.8 


9.39 29.7 


1.95 


2.61 


0.75 




0.67 0.24 


n 


-1.65 


-1.483 


17.1 


8.95 32.5 


0.84 


2.62 


0.55 




0.26 0.72 




-1.82 


-1.381 


22.3 


11.5 43.1 


0.73 


2.73 


0.00 




0.99 


4 


-1.44 


-1.327 


15.9 


8.30 30.4 


0.70 


2.59 


0.35 




0.96 



present result. We thus expect that the momentum distribution data of 15 C fragments from 
16 C breakup 16] are well reproduced in the present model as in the previous paper jij. 
The two-neutron correlation is examined with the density distribution function 

Pfru X2, 0) = 7pr~r E(^m(1, 2) | * JM (1, 2)) spin , (8) 

where 9 is the angle between x\ and x 2 and (• • -^pm indicates that the integration is to be 
done over the spin coordinates only. Figure |2] displays the contour map of 87r 2 x 4 sin6 l p(x, x, 6) 
for the + states. For the Of state, we have two distinct peaks: One peak with smaller angles 
is the highest (about 0.02 fm~ 2 in height), suggesting the correlation of "di-neutron" type, 
and the other with larger angles reaches half the highest peak, corresponding to a "cigar" - 
like configuration in which the two neutrons sit on the opposite sides of the core. For 
the 0^" state, however, only one distinct peak with a height of about 0.014 fm -2 appears 
around 6 = 90° (a "boomerang" shape). These different characteristics are expected from 




FIG. 2: Contour maps of the two-neutron densities p(x,x,9), weighted by 87r 2 x 4 sin#, for the 0^ 
and O2" states of 16 C. The contour shows peaks, and the difference between any two neighboring 
contour levels is 0.002 fm -2 . See text. 
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the expectation value of cos#, which is estimated from (x\'X%)/ (xf). This value is 0.17 for 
0^ and 0.06 for 0^, as seen from Table ITT1 We also confirmed that the 2f state has two 
peaks similarly to the ground state but the other states all have one peak at the boomerang 
configuration. 

The E2 operator, M.^ for 16 C in the 14 C+n+n model is expressed as 0] 

M„ = -Sey 2tl (r)+q eS ey 2tl (R), (9) 

with q e g— Jj + of^) where 5 is a polarization charge of the valence neutron. The value of 5 
can be estimated so as to fit the B(E2) value of 15 C within the 14 C+n model, which leads 
to 5=0.160 (set A), 0.159 (set B) and 0.151 (set C), respectively. The first term on the 
right-hand side of Eq. Q is the E2 operator corresponding to the relative motion of the two 
neutrons, while the second term corresponds to the relative motion between the 14 C core 
and the center-of-mass of the two neutrons. If the valence neutrons have no charge (5 = 0), 
the E2 operator reduces to the second term of Eq. (0) with a small qWj. In this case, the 
B(E2) value is contributed by the difference between the center of mass of 14 C and that of 
16 C. For nonzero 5, both terms of Eq. Q contribute to the E2 transition matrix element. 
The B(E2; 2f — > Of) value calculated with the different potentials are shown in Fig. El All 
the cases give similar results. The B(E2) value is about twice the measurement for such 5 
that is determined from the 15 C data. It is interesting to see how the B(E2) value changes 
with the admixture of the S = 1 components. The E2 operator (jUJ) is independent of spin, 
so its reduced matrix element consists of two parts: 

(0+||M||2+> = ^/(l- fl 2)(l-62) ( M (S = 0)) + ab(M(S = 1)), (10) 

where a and b are the probability amplitudes of the S = 1 components in the Of and 2f 
states, respectively, and the (A4(S)) stands for the reduced matrix element between their 
(normalized) wave function components with spin S. For set B and for 5 = 0.159, we have 
(M{S = 0)) = 2.75efm 2 and (M{S = 1)) = 1.56 efm 2 . If both of the Of and 2+ states are 
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5 

FIG. 3: The B(E2) value for the transition from the 2f state to the ground state in 16 C as a 
function of the polarization charge 5. The measured value P, 0] 0.63 ± 0.12 e 2 fm 4 is indicated by 
horizontal lines. 
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purely spin-singlet (a = b = 0), the B(E2) value would be 2.75 2 /5 = 1.51e 2 fm 4 . In fact, 
they have S = 1 admixtures whose magnitudes are a 2 = 0.08 and b 2 = 0.25 (see Table ITTj) . 
and therefore the B(E2) value reduces to 1.25 e 2 fm 4 , as shown in Fig. El In the previous 
calculation with S = only jj, {M{S = 0)) was 2.63 efm 2 , leading to the B(E2) value of 
1.38 e 2 fm 4 . Thus the S=l admixture has led to reducing the B(E2) value. 

We have studied the structure and anomalous E2 transition of 16 C in the 14 C+n+n model, 
where the two valence neutrons interact via the realistic potential. We have tested three 
different n- 14 C potentials which all reproduce the energies of the | x and | x states of 15 C in 
the 14 C+n model. Both spin singlet and triplet components are included in the calculation. 
The stochastic variational method has been used to solve the three-body problem with the 
Pauli constraint on the valence neutrons. The energy spectrum of 16 C below the 14 C+n+n 
threshold is in fair agreement with experiment except that the ground-state energy is about 
1 MeV too high. The dependence of the B(E2;2\ ^ 0^) value on the polarization charge 
of the valence neutrons is studied and found not to be very sensitive to the choice of the 
n- 14 C potential. We have seen that admixing the S = 1 components in the wave functions 
reduces the B(E2) value considerably. The calculated B(E2) value is, however, about twice 
the observed value if the polarization charge is set to reproduce the B[E2; | — ► | ) of 15 C. 

Within the 14 C+n+n model, the present calculation is probably the most unambiguous 
because it is free from any adjustable parameters and the result is rather insensitive to the 
choice of the n- 14 C potential. The fact that the ground-state energy is predicted to be 1 
MeV high and that the theoretical B(E2) value is still twice large may point to the necessity 
of other mechanisms such as core distortion or excitation. 
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